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Q^ , Abstract. In this invited talk I show how the reddest and rarest galaxies 

at high redshift {z ~ 1.5) can be used to set the stellar evolution clock. 
C , I argue that one can confidently compute the collapse redshift of these 

objects. This yields to a high collapse redshift [z > 6) and therefore 

their age is well constrained (in all cosmologies) between 3 and 4 Gyr. I 

also show that this is, indeed, the age derived using a variety of synthetic 

stellar population models when proper statistical tools are used to analyse 

^ ' their observed spectral energy distribution. This allows me to conclude 

^^ . that all stellar population models yield to the same consistent age for 

L-^ I these galaxies, i.e. about 3.5 Gyr and that the stellar clock is properly 

\^ • set. Low ages are therefore excluded with high confidence. 

o: 
o\. 

r^ ■ 1. Introduction 

Oh; 

Q ■ Traditionally, the stellar clock has been set on stellar objects at z = 0, namely the 

H \ Sun. Here I describe how one can find passively evolving galaxies at high-redshift 

t/^ ■ and use them to set the stellar clock of stellar populations of intermediate age. 
. . \ The most problematic issue is how to find galaxies at high redshift. The 

^^ ■ study of 'normal' star-forming galaxies at z > 2 has developed into a booming 

K^ , astronomical industry over the past 4 years (e.g. Steidel et al. 1998). Since 

^ • most high-redshift galaxies are optically selected, one is biased towards blue 
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objects that show recent star formation, i.e. biased towards composite stellar 
populations with several episodes of star formation. In this way, one can only 
hope to determine the age of the youngest stars at high-redshift, not a very 
useful age indicator. Nevertheless, a few valiant efforts have been carried out 
in order to determine the age of the stellar populations in high-redshift galaxies 
(Chambers & Chariot 1990). 

A more promising way to find passively evolving objects is utilizing radio 
galaxies. One of the cleanest results in extra-galactic astronomy is that all 
powerful (-P > 10^^ WHz'^sr"^) radio sources in the present-day universe are 
hosted by giant ellipticals. It is then reasonable to assume that high-redshift 
radio sources also reside in ellipticals or their progenitors. By selecting radio 
galaxies at mJy fiux levels we (Dunlop et al. 1996, Spinrad et al. 1997, Dunlop 
1998, Dey et al. 1999) have shown that that it is possible to find examples of well 
evolved galaxies at z ~ 1.5 whose near- ultraviolet spectrum is uncontaminated 
by a recent burst of star formation. Keck spectroscopy of these objects has 



yielded the first detection of stellar absorption features from old stars at z < 1.5 
and thus the first reliable age-dating of high-redshift objects. The best example 
in our sample is 53W069. Using WPFC2 and NICMOS images below and above 
the 4000 A break we have verified that the above also holds for its morphological 
properties and scalelengths (Dunlop 1998). Using a 2-dimensional fitting code 
it is possible to show that 53W069 is consistent with a r^'^ law and inconsistent 
with a exponential disc profile. Furthermore, a physical half-light radius of 
Te ~ 4 kpc has been obtained assuming Q = 1 and Hq = 50 km s^^ Mpc^^, 
which lies exactly in the Kormendy relation for ellipticals (Dunlop 1998). 

The traditional approach would be to use synthetic stellar population mod- 
els to derive the age of 53W069 from its observed spectral energy distribution. 
Here I argue, that it is possible to use observations of the abundances and cluster- 
ing of high-redshift galaxies to estimate the power spectrum on small scales, and 
thus constrain the age of 53W069 without the need of stellar evolution physics 
knowledge. The following section summarizes the results of this exercise, as 
given by Peacock et al. (1998). 

2. Small-scale po^ver spectrum 

2.1. Press- Schechter apparatus 

The standard framework for interpreting the abundances of high-redshift ob- 
jects in terms of structure-formation models, was outlined by Efstathiou & Rees 
(1988). The formalism of Press & Schechter (1974) gives a way of calculating 
the fraction Fc of the mass in the universe which has collapsed into objects more 
massive than some limit M: 



Fc{> M,z) = 1-erf 



V2a{M) 



(1) 



Here, a{M) is the rms fractional density contrast obtained by filtering the linear- 
theory density field on the required scale. In practice, this filtering is usually 
performed with a spherical 'top hat' filter of radius R, with a corresponding 
mass of AirpiyR^ /?>, where pi, is the background density. The number 5c is the 
linear-theory critical overdensity, which for a 'top-hat' overdensity undergoing 
spherical collapse is 1.686 - virtually independent of il. This form describes 
numerical simulations very well (see e.g. Ma & Bertschinger 1994). The main 
assumption is that the density field obeys Gaussian statistics, which is true in 
most inflationary models. Given some estimate of Fc, the number (t{R) can 
then be inferred. Note that for rare objects this is a pleasingly robust process: 
a large error in Fc will give only a small error in cr{R), because the abundance 
is exponentially sensitive to a. 

Total masses are of course ill-defined, and a better quantity to use is the 
velocity dispersion. Virial equilibrium for a halo of mass M and proper radius 
r demands a circular orbital velocity of V^ = ^^. For a spherically collapsed 
object this velocity can be converted directly into a Lagrangian comoving radius 
which contains the mass of the object within the virialization radius (e.g. White, 



Efstathiou & Frenk 1993): 






R/h-'Mvc=-^ — -^. (2) 



Here, Zc is the redshift of virialization; ilm, is the present value of the mat- 
ter density parameter; /c is the density contrast at virialization of the newly- 
collapsed object relative to the background, which is adequately approximated 
by /c = 178/ri5^^(zc), with only a slight sensitivity to whether A is non-zero 
(Eke, Cole & Frenk 1996). For isothermal-sphere haloes, the velocity dispersion 
is a^ = Vc/vl. Given a formation redshift of interest, and a velocity dispersion, 
there is then a direct route to the Lagrangian radius from which the proto-object 
collapsed. 

2.2. Abundances and masses of high-redshift objects 

Three classes of high-redshift object can be used to set constraints on the small- 
scale power spectrum at high redshift: 

(1) Damped Lyman-a absorbers 

If the fraction of baryons in the virialized dark matter halos equals the 
global value ri^B, then data on these systems can be used to infer the total 
fraction of matter that has collapsed into bound structures at high redshifts (see 
Peacock et al 1998 and refs. therein). Therefore, 

F, = ^ ~ 0.12/1 (3) 

for these systems. In this case alone, an explicit value of h is required in order 
to obtain the collapsed fraction; h = 0.65 is assumed and we have adopted 
nuBh'^ = 0.02. 

(2) Lyman- limit galaxies 

Steidel et al. (1996) identified star-forming galaxies between z = 3 and 3.5 
by looking for objects with a spectral break redwards of the U band. Steidel et 
al. give the comoving density of their galaxies as 

N{n = 1) ~ lO"^-^'^ ( h-^ Mpc)~^ (4) 

Direct dynamical determinations of these masses are still lacking in most 
cases. Steidel et al. attempt to infer a velocity width by looking at the equivalent 
width of the C and Si absorption lines. These are saturated lines, and so the 
equivalent width is sensitive to the velocity dispersion; values in the range 

o-„ ~ 180-320kms"^ (5) 

are implied. In practice, this uncertainty in the velocity does not produce an 
important uncertainty in the conclusions. 

(3) Red radio galaxies 

This is the set of observations for which we wish to determine their collapse 
redshift. Two extremely red galaxies were found at z = 1.43 and 1.55, over an 



area 1.68 x 10 ^ sr, so a minimal comoving density is from one galaxy in this 
redshift range: 

N{n = l)W~^-^^ (h-^Mpc)-^. (6) 

This figure is comparable to the density of the richest Abell clusters, and is 
thus in reasonable agreement with the discovery that rich high-redshift clus- 
ters appear to contain radio-quiet examples of similarly red galaxies (Dickinson 
1995). 

Since the velocity dispersions of these galaxies are not observed, they must 
be inferred indirectly. This is possible because of the known present-day Faber- 
Jackson relation for ellipticals. Their large- aperture absolute magnitude is 
Mv{z = 1.55 I 17 = 1) ~ —21.62 — 51ogio h (measured direct in the rest frame). 
This yields to a^ = 222 to 292 kms~^, which is a very reasonable range for a 
giant elliptical, and is adopted in the following analysis. 

Having established an abundance and an equivalent circular velocity for 
these galaxies, the treatment of them will differ in one critical way from the 
Lyman-a and Lyman-limit galaxies. For these, the normal Press-Schechter ap- 
proach assumes the systems under study to be newly born. For the Lyman-a and 
Lyman-limit galaxies, this may not be a bad approximation, since they are evolv- 
ing rapidly and/or display high levels of star-formation activity. For the radio 
galaxies, conversely, their inactivity suggests that they may have existed as dis- 
crete systems at redshifts much higher than z ~ 1.5. The strategy will therefore 
be to apply the Press-Schechter machinery at some unknown formation redshift, 
and see what range of redshift gives a consistent degree of inhomogeneity. 

2.3. Collapse redshifts and ages for red radio galaxies 

Fig. 1 shows the cr{R) data which result from the Press-Schechter analysis, for 
three cosmologies. The cr{R) numbers measured at various high redshifts have 
been translated to z = using the appropriate linear growth law for density 
perturbations. 

The open symbols give the results for the Lyman-limit (largest R) and 
Lyman-a (smallest R) systems. The approximately horizontal error bars show 
the effect of the quoted range of velocity dispersions for a fixed abundance; 
the vertical errors show the effect of changing the abundance by a factor 2 
at fixed velocity dispersion. The locus implied by the red radio galaxies sits 
in between. The different points show the effects of varying collapse redshift: 
Zc = 2, 4, . . . , 12 [lowest redshift gives lowest cr{R)]. Clearly, collapse redshifts of 
6-8 are favoured for consistency with the other data on high-redshift galaxies, 
independent of theoretical preconceptions and independent of the age of these 
galaxies{see also Kashlinsky &: Jimenez (1997)). 

What is then the age of the red radio galaxies as inferred by their high 
collapse redshifts? First bear in mind that in a hierarchy some of the stars in 
a galaxy will inevitably form in sub-units before the epoch of collapse. At the 
time of final collapse, the typical stellar age will be some fraction a of the age of 
the universe at that time: age = t{zohs) — t{zc) + at{zc). We can rule out a = 1 
(i.e. all stars forming in small subunits just after the big bang). For present- 
day ellipticals, the tight colour- magnitude relation only allows an approximate 
doubling of the mass through mergers since the termination of star formation 
(Bower at al. 1992). This corresponds to a ~ 0.3 (Peacock 1991). A non-zero a 
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Figure 1. The present-day linear fractional rnis fluctuation in den- 
sity averaged in spheres of radius R. The data points are Lyman-a 
galaxies (open cross) and Lyman-limit galaxies (open circles) The di- 
agonal band with solid points shows red radio galaxies with assumed 
collapse redshifts 2, 4, ... 12. The vertical error bars show the effect of 
a change in abundance by a factor 2. The horizontal errors correspond 
to different choices for the circular velocities of the dark-matter haloes 
that host the galaxies. The shaded region at large R gives the results 
inferred from galaxy clustering. The lines show CDM and MDM pre- 
dictions, with a large-scale normalization of cxg = 0.55 for 17 = 1 or 
(Tg = 1 for the low-density models. 
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Figure 2. Age of the Universe at z = 1.5 for several cosmologies, it 
transpires from the figure that the age of galaxies formed at z > 5 is 
well constrained between 3 and 4 Gyr 

just corresponds to scaling the collapse redshift as apparent (1+Zc) oc (1— q)^^'^, 
since t oc (1 + z)~^' ^ at high redshifts for all cosmologies. For example, a galaxy 
which collapsed at z = 6 would have an apparent age corresponding to a collapse 
redshift of 7.9 for a = 0.3. 

Converting the ages for the galaxies to an apparent collapse redshift depends 
on the cosmological model, but particularly on Hq. Some of this uncertainty may 
be circumvented by fixing the age of the universe. After all, it is of no interest 
to ask about formation redshifts in a model with e.g. = 1, /i = 0.7 when 
the whole universe then has an age of only 9.5 Gyr. If fi = 1 is to be tenable 
then either h < 0.5 against all the evidence or there must be an error in the 
stellar evolution timescale. If the stellar timescales are wrong by a fixed factor, 
then these two possibilities are degenerate. It therefore makes sense to measure 
galaxy ages only in units of the age of the universe - or, equivalently, to choose 
freely an apparent Hubble constant which gives the universe an age comparable 
to that inferred for globular clusters. In this spirit, Fig. 2 gives apparent ages 
as a function of effective collapse redshift for models in which the age of the 
universe is forced to be 14 Gyr (e.g. Jimenez et al. 1996, see also Fusi-Pecci in 
this volume). 

This plot shows that the ages of the red radio galaxies are not permitted 
very much freedom. Formation redshifts in the range 6 to 8 predict an age of 
close to 3.1 Gyr for i7 = 1, or 3.7 Gyr for low-density models, irrespective of 
whether A is nonzero. The age-Zc relation is rather flat, and this gives a robust 
estimate of age once we have some idea of Zc through the abundance arguments. 



3. The age of 53vi^069 

As seen from the previous section, galaxies at z ~ 1.5 with low comoving den- 
sities, are excellent sites for setting the stellar clock since their age is known 
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Figure 3. Left panel: the spectrum of 53W069 overlaid with the best 
fitting Bruzual-Charlot model, which has an age of 3.3 Gyr. Right 
panel: Same as before but now overlaid with the best fitting model 
from Jimenez et al. (1999b), which has an age of 4.0 Gyr. 



from other independent arguments. In what remains I will analyse how some of 
the different set of synthetic stellar population models available perform when 
trying to recover the age of 53W069 (the reddest passively evolving galaxy found 
at high— z). 

The spectral energy distribution (SED) of 53W069 is presented in Fig. 3. 
In order to determine its age I have performed properly weighted chi-squared 
fits to the ultra-violet SED using the popular Worthey and also Bruzual & 
Chariot synthetic stellar population models and the models developed by our 
group (Jimenez et al. 1999b). The results are listed in Table 1. A few important 
points transpire from this table. First, all models yield ages larger than 3 Gyr for 
53W069. Second, column 1 and 2 show that some models seem to be internally 
inconsistent, in the sense that they are capable of reproducing very red R — K 
colours at a much younger age than they can reproduce the ultraviolet SED or 
the spectral breaks. However, since R — K '\s mainly affected by the evolution of 
the late stages of stellar evolution, once should focus on ages derived from the 
UV-spectrum since this only depends on the correct prediction of the MSTO, 
a simpler and thus an easier part to model. Indeed, not do so it tantamount 
to throwing away the new, more robust information which can be gleaned from 
the spectroscopy. If one focuses on the results of 53W069 then, ignoring the 
anomalously young R — K age produced by some models, all sets are basically 
in good agreement that the overall shape of the UV SED is consistent with an 
age in the range 3.0 to 4.0 Gyr, and certainly yielding a robust (99% confidence) 
minimum age of 3.0 Gyr. Thus one finds that all models yield basically to the 
same result, in perfect agreement with the independently derived age in section 
2. Therefore, one can conclude that the stellar clock is properly set, even at these 
young ages. The anoumalously small ages (about 1 Gyr) obtained by fitting the 
R — K colours, can be entirely attributed to the fact that K (rest-frame 1 /u) is 



entirely dominated by the giants population, and thus by mass loss and other 
complications from stellar evolution during its late stages. 

Spectroscopically, 53W069 thus appears to be the best known example of 
old, passively-evolving elliptical galaxies at redshifts as high as z ~ 1.5. It is 
worth noting that z ~ 1.5 is the redshift where one expects galaxies to be redder, 
and not at higher redshift (see Jimenez et al. 1999a) 



Table 1. A comparison of age estimates for the stellar population of 
53W069 as derived from the instantaneous burst models of Bruzual & 
Chariot (B&C), Worthey (1998) (W) and Jimenez et al (1999b) (J99), 
when used to fit different spectral indicators of age. 



Feature 


B&C 


W98 


J99 


UV-SED 


3.3 Gyr 


3.1 Gyr 


3.7 Gyr 


R-K 


1.6 Gyr 


2.0 Gyr 


4.0 Gyr 


2649 A 


5.0 Gyr 


4.0 Gyr 


4.0 Gyr 


2900 A 


4.5 Gyr 


4.2 Gyr 


4.0 Gyr 



In summary, the new data on 53W069 clearly show that the Universe at 
z ~ 1.5 contains stellar systems whose populations are 3 to 4 Gyr old. At 
z ~ 1.5 the Universe was less than 30% of its present age, and the uncertainties 
are largely independent of those encountered in GCs studies. The existence 
of 53W069 permits only low Hubble constants and/or low cosmic densities; in 
particular, an $7 = 1 Universe requires ffo ^ 45 km s^^ Mpc^^ . On the other 
hand, in an Universe with a cosmological constant, 53W069 requires ^\ > 0.4 
if Ho > 60 km s"^ Mpc"^ 
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